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1. Introduction 

Recently high emphasis is placed on testing the unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix. In particular, the unitarity of the CKM matrix implies the unitarity 
constraint of the first row 

\V ud ? + \V us \ 2 +\V ub \ 2 = \-8. (1.1) 

According to the PDG [1] 8 is equal to 0.0008(11). The uncertainty is still quite substantial and 
has to be decreased. As \V u b\ is much less than unity, about half of the error of 8 comes from 
the uncertainty of \V US \. Since the kaon semileptonic decay rate is proportional to |V iU | 2 |/ + (0)| 2 , 
this matrix element can be determined by combining experimental results for this decay rate and 
theoretical calculations of the vector form factor at zero momentum transfer /+(0). 
The kaon semileptonic decay form factors f+(q 2 ) and f-{q 2 ) are defined as 

{K{p')\V^\K{p))=f + {q 2 ){p + p') ll +f^ q 2 ){p-p')^ (1.2) 

where q 2 = (p — p') 2 is the momentum transfer and Vn = sy^u is the weak current. It is also 
convenient to introduce the scalar form factor defined as 

A <« 2 > = (' + 1^^) ■ ™ d ^> - tM ■ <L3) 

The form factor /+(0) at zero momentum transfer was estimated by Leutwyler and Roos [2] in 
1984. They obtained the value 0.961(8), which is still used as a reference value to extract \ V US \ from 
the experimental data. However, to estimate higher order terms in the chiral perturbation theory 
(ChPT) expansion, they used a model of the wave function of the pseudoscalar meson. Recent 
ChPT calculations [3] favour a slightly larger value of /+(0) = 0.984(12). Thus, it is desirable to 
calculate /+ (0) non-perturbatively and lattice calculation provides such an opportunity. 

Recently a lot of lattice groups reported calculations of /+(0) [4, 5, 6, 7, 8, 9, 10, 11]. Their 
results coincide with each other and are in agreement with that of Leutwyler and Roos. For a recent 
review see [12]. 

In this paper we present first preliminary results for these form factors from Nf = 2 flavours 
of light dynamical, non-perturbatively 0(a)-improved Wilson fermions and Wilson glue. 

2. Lattice setup 

For our calculation we use a gauge ensemble at /3 = 5.29, which conesponds to a lattice spac- 
ing a = 0.075 fm, 1 and fc sea = 0.13590, which corresponds to a pion mass M K = 0.591 (2) GeV. 
The lattice volume is 24 x 48 with spacial extent equal to 1.9fm. The correlation functions have 
been calculated on about 800 gauge configurations using various source locations to reduce statis- 
tical noise. We identify the sea quark mass with the light quark mass and the valence quark mass 
with the strange quark mass using 3 values K s = 0. 13485, 0. 13530 and 0. 13570. The corresponding 
kaon masses are M K = 0.780( 11), 0.704( 13) , 0.629( 14) GeV. 

'For translation into physical units the Sommer parameter tq = 0.467 fm is used. 
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3. Correlators 

On the lattice we calculate two and three point functions, which are defined by 

C P (t,p) UO P . snk (xj)ol rc (0,0))e-^ > Z ^yV W (3.1) 

Y ' >^°° 2Ep(p) 

Z * p >° Asnk (e(p')|vi* ) |P(p))e- & ^ , -^»') ( *'- , >, (3.3) 



f,(t'-0— 4E P {p)E Q {p')Z v 

where P and 2 denote either the or tt meson. The energy of meson P (Q) is denoted by Ep(p) 
(En(p)). The renormalised vector current, including the renormalisation factor Zy, is denoted by 

(R) 

Vn ■ The overlap with the physical meson states is given by Zp^ and Z PiSrc . 

In the following we will assume the point meson sources to be inserted at time t src = (thus 
t src is omitted in Eq. (3.3)) and point sinks are inserted at t' = T/2. As t' is fixed the t' dependence 
of all quantities is ignored in the following. For this choice of t' the three point functions (and 
therefore also the double ratios) are symmetric with respect to T/2. We make use of this property 
and average over both time ranges to increase the precision of calculations. 

Note that smearing of meson operators is not used. This is because it leads to a momentum 
dependence of the overlap Zp^ sn k{p) and Zp sn {p). See, for example, Appendix A in Ref. [13] and 
the discussion below. 



4. Scalar form factor at q 2 max 

The scalar form factor fo(q ) at q 2 = qj nax = (Mg — M K ) 2 can be obtained from the double 
ratio of three point functions (which was originally proposed to calculate the B — ► Dlv form factor 
in Ref. [14]): 

_ C««(t,0,d)Cf(tm {M K + M n ) 2 2 2 

R{t) -CfHtA0)C^ t ,0,0)^ 4M A (/0(w)) -*■ 

Our results for R(t) for three values of the strange quark masses are shown in the left panel of 
Fig. 1. Note that renormalisation constants and exponential factors exactly cancel in R(t). We can 
extract fo(q 2 nax ) with a statistical uncertainty < 0.1%. 

In the SU (3) symmetric limit R is equal to unity. The deviation of R from unity depends on 
the physical SU (3) breaking effects on /ot), which can be seen in the right panel of Fig. 1. 

5. Interpolation to zero momentum transfer 

To study the q 2 dependence of the form factor we calculate 

_ ^, f+{q 2 ) ( E K {p)-E n {p') 2 \ 
F{p,p)= f , 2 , 1+ ,^ - ,~,M q)) (5-1) 

Mqlax) V E K (p)+E n {p') J 



3 



Kaon semileptonic decay form factors 



S.M. Morozov 



1.04 



1.03 



-1 .02 ; 



k s = 0.13485 
k_ = 0.13530 
k q = 0.13570 



1.01 - } 5 5 J 



i i i i i i 



i i i i 



1.00 



10 , /a 15 20 



1 .0025 
1.002 
-a 1.0015 

rd 
E 

CM 

-° 1.001 
1 .0005 
1 



0.01 



0.02 0.03 

a 2 (M K 2 - M„ 2 ) 



0.04 0.05 



Figure 1: Left: Time dependence of the double ratio R(t) (Eq. (4.1)) for three values of the strange quark 
masses. Right: Values of fo(ql wx ) as a function of the SU (3) breaking parameter a 2 AM 2 = a 2 (M^—M 2 .). 

from the second double ratio 

K(tit ^ Cf*(t,p,p')C K (t,d)C*(T/2-t,0) E K (p)+E n (p>) 

Kp(t,p,p ) = ^-=i > r(P,P ) = Rp . (5.2) 

1 ' Cf*(t,0,0)C K (t,p)C*(T/2-t,p>) M K + M n ' 

For meson point operators the overlaps cancel in (5.2). This is not true anymore if smeared 
meson operators are used, because their overlap with the physical state depends on the momentum. 
That is why we are using point sources and sinks. A systematic study of this effect will be reported 
elsewhere. 

The double ratio Rp{p,p') is symmetric under exchange of 71 and K and their momenta, i.e. 
Rp~^ n {p,p') = Rp^ K (p' ,p), which we make use of to increase the statistics at some values of q 2 . 
In the left panel of Fig. 2 we show, as an example, the data for Rp(t,p,p') with \p'\ = 0, \p\ = 1 
for K —>■ % and \p'\ = 1, \p\ = for the % — > K transition. 

In order to convert F(p,p') to f + (q 2 ), we need to calculate t, (q 2 ). To do so, we define a third 
double ratio: 



^ C**(t,p,p')CF(t,p,p') 



(5.3) 



Then ^ (q ) is given by 



h 2 ) = -(p+p')k(E K (p)+E K (p')) + (p+p') k (E K (p)+E4p'))R k (p,p') 54 

SW ' {p-p%{E K {p)+E K {p'))-{p + p') k {E K {p)-E 7l {p'))R k {p,p') ' 1 ■ 

where R k (p,p') = limf-,™ Rk(t,p,p'). The double ratio R k (p,p') is exactly unity in the SU(3) 
symmetric limit and is sensitive to SU (3) breaking effects in the matrix element Cjp 5 . 

This ratio is the noisiest of all of them. Thus, ^(q 2 ) could only be extracted with a large 
uncertainty. The relative error turns out to be 20% — 80% for our statistics, even after we make use 
of the symmetry: 

Rf-(p,p')=Rr K (p',p)=Rr K (-p',-p). 
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Figure 2: Left: Time dependence of double ratio Rp(t ,p,p ! ) for K s — 0.13485. Right: Scalar form factor 
fo(q 2 ) for different values of the strange quark mass. The solid line is the result of a monopole fit as described 
in the text. 

From fo{q 2 nax ), F(p,p') and t, (q 2 ) we can calculate the scalar form factor fo(q 2 ). The results 
are shown in the right panel of Fig. 2. To interpolate fo(q 2 ) to zero momentum transfer, we fit 
it with a monopole ansatz fo(q 2 ) = fo{0) /(l — q 2 /M 2 ). The result of the fit is also shown in the 
figure. In the given range of momenta the data is well described by this ansatz. 

6. Chiral extrapolation 

In order to calculate the physical value of /+(0) = /o(0) we have to extrapolate our results to 
the physical pion and kaon masses. We make use of the results of ChPT to guide our extrapolation. 
In ChPT /+(0) can be expanded in terms of light pseudoscalar meson masses giving 

/+(0) = l+/ 2 +/4 + ... ,fn = 0(Ml n K>J1 ). (6.1) 

The leading correction f% receives only contributions from non-local operators and can be deter- 
mined unambiguously in terms of Mk,M k and f n (see [15]). We compute fi at the actual pion and 
kaon masses and define 

A/ = /+(0)-(l+/ 2 ), (6.2) 

which receives only contributions from local operators. 

The Ademollo-Gatto theorem [16] states that Af is proportional to {M\ —M 2 ) 2 . Hence, we 
may write 

Af = a + b(M 2 -M 2 n ) 2 , (6.3) 

to obtain /+(0) by extrapolating (6.3) to the physical point. The data points and the resulting fit are 
shown in the left panel of Fig. 3. We find a « 0, in agreement with the expected result Af = in 
the limit of flavour SU (3). 
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Figure 3: Interpolation of A/ (Left) and | (0) (Right) to physical meson masses indicated by the vertical 
lines. 

At the physical meson masses we find Af = —0.0126(15) 2 . This is to be compared with 
Af = -0.016(8) of Ref. [2]. Inserting the physical value of f 2 , f 2 = -0.0227, into Af, we obtain 

/+(0) = 0.9647(15 ) ttat . (6.4) 

We extract the CKM matrix element \V US \ from the experimental value of |V Mi |/+ (0) averaged over 
all decay modes [17], \V us \f+(0) = 0.21673(46). This finally gives 

\V US \ =0.2247(5)^(4)^,. (6.5) 

To compute <^(0) at the physical meson masses, we make the ansatz ^(0) = c{M\ — M%), in 
accord with t, (0) = in the SU (3) limit. In the right panel of Fig. 3 we show the data points and 
the resulting fit. We obtain 

£(0) = -0.10(2)^. 

This value is consistent with the experimental values —0.01(6) from Kf 3 decay and —0.125(23) 
from decay. 

Note that in our analysis only statistical errors are estimated. Analysis of systematic errors 
which include extrapolation errors of fo{q 2 ) to zero momentum transfer and chiral extrapolation 
errors is underway. 

7. Conclusions and outlook 

In this work we have presented preliminary results for the kaon semileptonic decay form fac- 
tors and \ V US \ from Nf = 2 non-perturbatively 0(a) -improved Wilson fermions. We found 

/+(0) = 0.9647(15)** , \V US \ = 0.2247(5)^(4)^ , (7.1) 

2 Here and in the following only the statistical error is quoted. 
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in agreement with the results of other lattice groups [4, 5, 6, 7, 8, 9, 10] as well as the estimate of 
Leutwyler and Roos [2], within the error bars. 

The QCDSF collaboration is going to improve the accuracy of the results by performing calcu- 
lations at lighter quark masses down to pion masses of about 300 MeV and using partially twisted 
boundary conditions [9]. In addition, we are going to study discretisation effects using gauge en- 
sembles at different lattice spacings in the range of 0.1 15 fm to 0.070 fm. Furthermore, we plan to 
study finite size effects. 
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